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ABSTRACT 

An analysis of data from the Spitzer Space Telescope, Hubble Space Tele- 
scope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is 
presented for the z = 0.036 merging galaxy system II Zw 096 (CGCG 448-020). 
Because II Zw 096 has an infrared luminosity of log(Lm/ L®) = 11.94, it is clas- 
sified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the 
Great Observatories All-sky LIRG Survey (GOALSj). The Spitzer data suggest 
that 80% of the total infrared luminosity comes from an extremely compact, red 
source not associated with the nuclei of the merging galaxies. The Spitzer mid- 
infrared spectra indicate no high-ionization lines from a buried active galactic 
nucleus in this source. The strong detection of the 3.3 /im and 6.2 jum PAH emis- 
sion features in the AKARI and Spitzer spectra also implies that the energy source 
of II Zw 096 is a starburst. Based on Spitzer infrared imaging and A KARI near- 
infrared spectroscopy, the star formation rate is estimated to be 120 M yr _1 and 
> 45M yr~ 1 , respectively. Finally, the high-resolution B, I, and .frband images 
show many star clusters in the interacting system. The colors of these clusters 
suggest at least two populations - one with an age of 1 — 5 Myr and one with 
an age of 20 — 500 Myr, reddened by — 2 magnitudes of visual extinction. The 
masses of these clusters span a range between 10 6 — 10 8 M & . This starburst source 
is reminiscent of the extra-nuclear starburst seen in NGC 4038/9 (the Antennae 
Galaxies) and Arp 299 but approximately an order of magnitude more luminous 
than the Antennae. The source is remarkable in that the off-nuclear infrared 
luminosity dominates the enitre system. 

Subject headings: infrared: galaxies — galaxies: individual(II Zw 096) — galaxies: 



^ttp:/ /goals. ipac.caltech.edu 
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interactions — galaxies: starburst — galaxies: star clusters 
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Introduction 



II Zw 096 (also known as CGCG 448-020 or IRAS 20550+1655) has an infrared 



lumin osity of log(Lj#/L ) = 11.94 and a luminosity distance of 161 Mpc ( lArmus et al. 
20091 ). Since its luminosity is above Lir > 10 11 L Q , it is classified as a Luminous Infrared 



Galaxy, or LIRG. From its optical morphology, II Zw 096 appears to be a merger of at least 



Arri 



two gas-rich spirals. From the optical imaging 
Class III interacting galaxy, based on the system of 
two identifiable nuclei, with well developed tidal tails 



)as et al. 



Suracd ( 



(120041) classify II Zw 096 as a 



19981 ). meaning the galaxy has 



With near-infrared imaging and spectroscopy, iGoldader et al.l ( 119971 ) uncovered two 
extremely red sources to the east of the merging disks that appeared to mark the location of 
a highly obscured, young starburst. These authors suggest that II Zw 096 is one of a handful 
of LIRGs (such as Arp 299 and VV 114) that are experiencing enhanced star- formation 



before the final dissipative collapse of the syste m. Whi 



observed spectroscopically in the mid-infrared ( IDudley 



e II Z w 096 has previously been 



19991 ). the ground based data 



were not of high enough signal-to-noise to permit detection of PAH or the fine structure 



lines. It is also known that II Zw 09 6 host s an O 
R.A.=20 h 55 m 05 s .3, Dec=+16°56'03" (IBaan 



1989). 



i megamaser at t he B 



Baan et al. 



(I1998h and 



950 position of 



Baan fc Klockner 



(120061 ) classified this megamaser as a starburst using the optical line ratios and the radio 



data. 

As part of the Great Observatories All-Sky LIRG Survey (GOALS), we have obtained 
images and spectra with the Spitzer Space Telescope, the Hubble Space Telescope (HST), 
the Chandra X-ray Observatory, and the Galaxy Evolution Explorer (GALEX) of a complete 
sample of LIRGs in the local universe. The G OALS targets consist of all LIRGs found in the 



IRAS Revised Bright Galaxy Sample (RBGS; 



Sanders et al. 



(120031 ) ). which covers galactic 



latitudes greater than five degrees and includes 629 extragalactic objects with 60 /zm flux 
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densities greater than 5.24 Jy. The median and the maximum redshift of the IRAS RBGS 
are z=0.008 and 0.088, respectively. The GOALS sample includes 179 LIRGs and 23 Ultra 
Luminous Infrared Galaxies (ULIRGs, Lm > 10 12 L Q ), covering the full range of galaxy 
interaction stages from isolated spirals to late stage mergers. GOALS provides an excellent 
dataset with which to explore the effect of mergers on infrared activity at low redshift. A 
critical part of the GOALS survey has been to use the Spitzer Space Telescope to identify 
the location of the infrared emission within LIR Gs, and to characte rize the source of the 
power. The GOALS project is fully described in 



Armus et al. 



((20091) 



Here, we present the first mid and far-infrared imaging and spectroscopy of II Zw 096 
from Spitzer, together with an analysis of the far-ultraviolet and optical imaging (HST), 
X-ray imaging (Chandra), and near-infrared spectra (AKART). In §2j we describe the 
observations and data reduction, in §21 we present our results, identifying the location of 
the bulk of the far-infrared emission, and in §H we discuss the nature of this buried power 
source. Cosmological parameters Hq = 70kms _1 Mpc -1 , fl m = 0.28, and Qa = 0.72 are 
used throughout this paper. 



Observations and Data Reduction 



2.1. Spitzer Imaging 



Spitzer observ ations of II Zw 096 were taken with the Infrared Array Camera (IRAC; 



Fazio et al. 



(12004 ) at 3.6 /^m, 4.5 /mi, 5.8 /mi, and 8.0 um on 2 04 O ctober 29 and the 



Rieke et al. 



feoOJ)) at 24 /mi, 70 /mi, 



Multiband Imaging Photometer for Spitzer (MIPS; 
and 160 /mi on 2004 November 30 (PID 3672, PI J. Mazzarella). All IRAC data were 
collected with 5 frames of 30 sec high dynamic range (HDR) mode. The 8.0 /mi data 



saturated the 30 sec integrations, so only the short 1.2 sec frames were used. Total exposure 
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times were 150 sec at 3.6 /im, 4.5 /xm, and 5.8 /xm and 6 sec at 8.0 /xm. The MIPS images 
were taken in photometry mode, using super-resolution at 70 /xm. Sixteen images of three 
seconds each were obtained. The total exposure times were 48.2 sec, 37.7 sec, and 25.2 sec 
for 24 /xm, 70 /xm, and 160 /xm, respectively. A summary of the Spitzer observations is given 
in Table [TJ 

We used the MOPExP software package to process the IRAC and the MIPS images, 
correcting the background, aligning, resampling the images, removing bad pixels, combining 
into mosaics, and clipping cosmic-ray events. Final mosaics were constructed with a pixel 
size of 0.6"/pixel for the IRAC images, 1.8"/pixel at 24 /xm, 4.0"/pixel at 70 /xm, and 
8.0"/pixel at 160 /xm for the MIPS images. 

The WCS of the images are limited to ~ 0.4" by the telescope pointing accuracy. 
However, alignment of the IRAC images has been refined with the Two Micron All Sky 
Survey (2MASS) to achieve a final accuracy (< 0.3"). For the MIPS images, the positional 
error is approximately 0.7" based on the startracker-to-boresight uncertainty and the MIPS 
scan mirror uncertainty (Lacy M. et al. 2009, in preparation). A pointing refinement with 
2MASS cannot be done because of the large wavelength difference between the MIPS and 
2MASS bands. 

All photometric fluxes were derived using the IDL APER routine. For the IRAC 
images, we u s ed 1.4 " and 2.9" radii circular apertures, using the aperture corrections from 



Surace et al. 



( 120041 ) . For the MIPS images, we fit the PSF to the peak of emission to 



estimate total flux. 



2 Available from the Spitzer Science Center 

http:/ /ssc. spitzer. caltech.edu/postbcd/mopex. html 



see 



-9 - 



2.2. HST 

The HST far-ultraviolet (FUV) and optical images were obtained with the Solar 
Blind Channel (SBC) and the Wide Field Camera (WFC) on the Advanced Camera for 
Surveys (ACS) on 2008 May 1 (PID 11196, PI A. Evans) and 2006 April 15 (PID 10592, 
PI A. Evans), respectively. The data were taken in the ACCUM mode using the LINE 
dither pattern for the WFC and the PARALLELOGRAM dither pattern for the SBC. The 
near-infrared images were taken using the NIC2 camera on the Near Infrared Camera and 
Multi-Object Spectrometer (NICMOS) with the MULTIACCUM mode and the SPIRAL 
dither pattern on 2008 June 1 (PID 11235, PI J. Surace). The following filters were 
used: ACS/SBC F140LP (FUV), ACS/WFC F435W (5-band), F814W (/-band), and 
NICMOS/NIC2 F160W (F-band) for 2528, 2520, 1440, and 2304 sec, respectively. We 
summarize the HST imaging observations in Table [TJ 

The HST data reduction began with the calibrated data products produced by STScI. 
The ACS data were reduced with the PyDrizzle software included in IRAF/STSDAS 
provided by STScI, to identify and reject cosmic rays and bad pixels, to remove geometric 
distortion, and to combine the images into mosaics. The NICMOS images were reduced 
using IRAF and IDL tasks designed to remove cosmic rays and bad pixels, and correct for 
the different bias levels across the quadrants due to the four amplifiers. The alignments 
of the HST images were improved based on the 2MASS astrometric solution. IRAF was 
used to construct the final mosaic, and the IDL routine APER was used to measure the 
photometry. An aperture size of 0.25", with corresponding aperture corrections of 1.18, 
1.16, 1.16, and 1.39, were used for photometric measurements of point sources in the final 
ACS/SBC F140LP, ACS/WFC F435W, F814W, and NICMOS F160W images, respectively. 



2.3. Chandra 



X-ray observations of II Zw 096 were obtained with the AXAF CCD Imaging 
Spectrometer (ACIS) on the Chandra X-ray Observatory with VFAINT mode on 2007 
September 10 (PID 8700551, PI D. Sanders). The background counts for total integration 
times 14.56 ks are 0.0025 counts/pixel and 0.004 counts/pixel in the soft (0.5 — 2 keV) 
and the hard (2 — 7 keV) bands, respectively. The 5-cr depths of the images are 
1 x 10~ 15 erg s _1 cm -2 in the soft band and 5 x 10 -15 erg s _1 cm -2 in the hard band. The 
fundamental data reduction was done by the Chandra data analysis software CIAO version 
3.4 and we have verified that there are no anomalous objects left in the data. Then we 
selected soft (0.5 — 2 keV) and hard (2 — 7 keV) X-ray energy ranges for making images and 
smoothed the data using a Gaussian filter with a a of 1.5 pixels (FWHM ~ 1.8"). 



2.4. Spitzer IRS Spectra 

II Zw 9 6 was observed with all modules of the Spitzer Infrared Spectrograph (IRS; 



Houck et al. 



(120041 )) in staring mode on 2006 November 12 (PID 30323, PI L. Armus). 
The slits were centered on the peak of the 8 /xm emission. The slit positions for all four 
modules are overlaid on the IRAC 8 jum image in Figure [H The number of nod cycles and 
the integration time in each nod position are listed in Table [TJ 

The data were reduced using the S15.3 IRS pipeline at the Spitzer Science Center. 
The pipeline software removes bad pixels and droop, subtracts the background, corrects 
linearity, and performs a wavelength and flux calibration. The backgrounds were subtracted 
with dedicated sky observations for the high-resolution data, and with off-source nods for 
the low-resolution data. One dimensional spectra were extracted with the SPICljfl software 



3 http: / / ssc.spitzer.caltech.edu/postbcd/spice.html 
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package, using the standard point source extraction apertures for all slits. For the high-res 
slits, this is equivalent to a full-slit extraction. The slit widths of Short-Low (SL), Long-Low 
(LL), Short-High (SH), and Long-High (LH) are 3.6", 10.5", 4.7", and 11.1", respectively. 
The extraction size for SL is 9.5", LL is 38.3", SH is 11.7", and LH is 22.8". 



2.5. AKARI Near-Infrared Spectra 

The near-infrared spectra of II Zw 096 were taken with AKARI on 2009 May 12 and 13 
(PID GOALS, PI H. Inami) durin g the post-helium (warm phase) mission using the Infrared 



Camera (IRC; 



Onaka et al. 



(120071 ) ) . The 2.5 — 5 /im spectra were obtained using the 1' x 1' 
aperture (the Np aperture) using a moderate-resolution grism. The spectral resolving 
power is ~ 120 at 3.5 /mi for a point-like source (~ 5" FWHM). The use of the Np aperture 
is intended to avoid the confusion of spectra but multiple sources in this aperture may 
interfere each other (see §3.5p . The observations used the IRCZ4 Astronomical Observation 
Template, taking four spectra, a reference image, and then five spectra at the end of the 
sequence. The exposure time for each spectrum was 44.41 seconds. Two pointings were 
used (in total 2 pointings x 9 frames) for the observations and the total exposure time was 
799 seconds. 

The data reduction was performed with the standard software package provided by 
JAXA, the IRC Spectroscopy Toolkit for Phase 3 data Version 20090211b This software 
performs dark subtraction, linearity correction, flat correction, background subtraction, 
wavelength calibration, spectral inclination correction and spectral response calibration. 
We extracted the ID spectra from the 2D spectral image using an aperture three pixels 
wide (4.4"). The number of bad pixels is increasing dramatically in the post-helium mission 



4 This is available from 



http://www.ir.isas.ac.jp/ASTRO-F/Observation/DataReduction/IRC/ 
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due to the elevated detector temperatures. A 3 x 3 moving window was used to isolate 
and replace bad pixels in each spectral image before co-adding the data and extracting the 
ID spectrum. The wavelength accuracy is evaluat ed to be 0.01 /im and the uncertainty in 



absolute flux calibration is estimated to be ~10% rtOhyama et al. 



20021). 



Results 



3.1. Spitzer Imaging 

In Figure [TJ the Spitzer images from IRAC 3.6 /im through MIPS 70 /im are presented 
in order of increasing wavelength. Three distinct objects are seen in the IRAC data. 
he nuclei of the me rging galaxies are to the north (source B) and south (source A) 



( iGoldader et al 



19971 ). In addition there is a bright, red source to the east (sources C and 
D), which is resolved and elongated at an angle of roughly 60°. This resolved IRAC source is 
coincident with the near-infrared sources "C" (R.A =20 h 57 m 24 s .47, Dec=+17°07'39 9") and 



Goldader et al 



fll997f ). We estimate the 



"D" (R.A.=20 h 57 m 24 s .34, Dec=+17°07'39 1") from 
flux density of source D at 3.6 /tm, 4.5 /im, 5.8 /im, and 8.0 /im, as 2.85 ± 0.08, 4.67 ± 0.13, 
12.7 ± 0.2, and 41.9 ± 0.8 mJy, respectively, using a circular aperture size of 1.4" radius. 

The MIPS 24 /im image shows a strong point source coincident with the position of 
source D, along with a lot of emission to the southwest, coincident with source A. Although 
the MIPS 70 jum image is unresolved, the position of the peak is consistent with source 
D. PSF fitting to source D at 24 /im and 70 /im suggests that it produces 67.2 ± 10.6% 
(1.37 ± 0.18 Jy) and 87.9 ± 13.8% (10.2 ± 1.2 Jy) of the total galaxy flux at 24 /im and 
70 /im, respectively. The angular separation of sources C and D is ~ 2.0" (see Figure [2] or 
H]). Although this small separation cannot be resolved with MIPS, the centroid of source D 
can be determined much more accurately than the FWHM of the beam. In addition, we 



13 



have simulated the MIPS image using two point sources separated by 2.0", determining that 
a flux ratio of 4 : 1 (D : C) is required to match the data. The photometry is summarized 
in Table gj 

We estimate the total infrared (TIR, 3 — 1100 /im) luminosity of source D to be 



TIR = 5.85 x 



11 L m from the Spitzer 24 and 70 /im images, using the formula Eq. (4) in 



Dale &: Heloul (120021 ) . assuming that the 160 /im flux of source D is 87.9% of the total at 



70 yum. This is approximately 80% of the TI R of the system. I I Zw 096 has an infrared 
luminosity (LIR, 8 - 1000 /im) of 1O 1194 L (jArmus et al.ll2009[ ). Assuming source D 
accounts for ~ 80% of the total LIR of the system (the same fraction as TIR), we estimate 



an LIR of 6.87 x 10 11 L Q . 



3.2. HST Imaging 

HST/ ACS and NICMOS images of II Zw 096 are shown in Figure [H from the upper 
left to the lower right, ACS/SBC F140LP (FUV), ACS/WFC F435W (5-band), F814W 
(/-band), and NICMOS/NIC2 F160W (//-band). The contours of the Spitzer 24 /im image 
are overlaid in each panel. 

The /-band and the //-band images with the 24 jum contour overlayed in Figure [2] 
suggest that source D is the origin of the bulk of the far-infrared emission in II Zw 096. 
The contour of the MIPS 70 /im image is also consistent with the position of source D. 

The HST NICMOS images resolve the structure of the red sources into at least 
10 "knots" (see Figure E]). Most of these knots are not seen in the visible image due 
to heavy extinction. In the //-band, the absolute magnitudes of sources C and D are 
Mjj = —18.7 ± 0.1 mag and —19.5 ± 0.1 mag, respectively, measured with an 0.25" radius 
aperture. We estimate that approximately 0.81 ± 0.04% and 1.62 ± 0.06% of the //-band 



-14- 



flux of the entire system comes from sources C and D alone. Sources C and D provide about 
7.78 ± 0.39% and 15.5 ± 0.6% of the if -band emission for the red complex in Figure EJ 
respectively. 

In Figure HJ we show a false-color RGB image of II Zw 096 made from the NICMOS 
if -band, ACS -B-band and ACS FUV images. The extremely red sources C and D are 
obvious in this figure, as are a number of unresolved blue sources to the west. These are 
presumably young star clusters in the northern spiral arm of the western galaxy and the 
overlap region. 

3.3. Chandra Imaging 

The Chandra X-ray images are shown in Figure [51 In the full (0.5 — 7 keV) and soft 
(0.5 — 2 keV) bands, source A dominates the emission. In the hard-band (2 — 7 keV) image, 
sources C and D become much more prominent. The 0.5 — 2 keV luminosity of source A 
and source C+D complex are 6.6 x 10 40 ergs _1 and 1.8 x 10 40 ergs _1 , respectively. In the 
hard band, the 2 — 7 keV luminosity of source A is 1 x 10 41 ergs -1 and the total of sources 
C and D is 1 x 10 41 ergs _1 . 

The C+D complex is clearly extended (see Figure [5]). Source C dominates the soft-band 
emission from the complex, with the fraction of the X-ray emission coming from source D 
apparently rising towards the hard band. This is evidenced by a gradual shifting of the 
peak from NE to SW in Figure [51 However, given the low number of total counts (~ 20) 
from this region and the strong overlap of source C and D, it is difficult to derive an 
independent measure of the hardness ratio of C and D, separately. Together, the hardness 
ratio (Hard-Soft/Hard+Soft) of source C+D is 0.05 ± 0.14. 
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3.4. Spitzer IRS Spectroscopy 

The Spitzer IRS Short-Low (SL), Long-Low (LL), Short-High (SH), and Long-High 
(LH) spectra of sources C and D (C+D) are shown in Figure [61 and the fluxes of the 
detected lines are shown in Table [2] and [3J Note that the LL and LH spectra cover not only 
source C+D but also source A due to their limited resolutions (see figured] (d)). However, 
the scale factor of LL/SL is 1.09, indicating the continuum is dominated by source C+D 
only. The IRS low-resolution spectra are dominated by strong PAH emission and silicate 
absorption. The high-resolution data show strong fine structure line emissions such as 
[Ne III], [Ne II], and [S III]. 

From the SL data, we measure 6.2 /un and 11.2 /im PAH equivalent widths (EQW) 
of 0.26 and 0.14 /im, respectively, and line fluxes of (5.57 ± 0.19) x 10~ 16 Wm~ 2 and 
(3.24 ± 0.32) x 10 _16 Wm" 2 , respectively, using a simple spline fit to the local continuum 
under the features. The apparent silicate optical depth at 9.7/xm is T9.7 ~ 1 when a screen 



geom etry is assumed. This implies an Ay > 19 mag toward source C+D (jRoche &: Aitken 



1984 ). This is com parable with the estim ate of 1 — 2 mag in the i^-band from the 



near-infrared colors (IGoldader et al 



19970 



We do not see evidence for 14.3 /im [Ne V] emission in the SH spectrum. A 3-er 
upper limit of the [Ne V] line flux is < 0.10 x 10~ 16 Wm~ 2 . The [Ne II] line has a flux of 
(2.74 + 0.47) x 10~ 16 W ur 2 , implying a 3-a upper limit on the [Ne V]/[Ne II] line flux ratio of 
0.04. The 25.9 /im [O IV] line is not detected, with a 3-a upper limit of < 0.29 x 10~ 16 W m" 2 , 
implying an upper limit on the [O IV] /[Ne II] line flux ratio of 0.11. The 18.71 /im [S III] 
line fluxes are detected in both of SH and LH: (2.16 ± 0.22) x 10~ 16 Wm" 2 in SH and 
(3.27 ± 0.04) x 10~ 16 Wm" 2 in LH. The large value for 18.71 /im [S III] in LH implies that 
there is measurable emission line flux from source A. The flux of the [S III] line at 33.48 /im 
is (4.77 ± 0.07) x HT 16 Wm" 2 . Using the [S III] 18.71/[S III] 33.48 line flux ratio in LH 
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of 0.69 ± 0.01, we derive an ionized gas density of 250 cm 3 at 10 4 K. The density is an 
average value over the area enclosed by LH which includes sources C, D, and A. 



-2 



The fluxes of H 2 S(3), H 2 S(2), and H 2 S(l) lines are (0.45 ± 0.05) x 10" 16 W m" 
(0.29 ± 0.03) x 10~ 16 W m" 2 , and (0.73 ± 0.05) x 10~ 16 W m" 2 , respectively. Taking these 
lines, and assuming an ortho to para ratio of 3.0, we estimate the warm molecular gas has 
a temperature of 329 K, and a total mass of 4.5 x 10 M & . 



f|2007h ). we 



If we fit the entire IRS low-res spectrum using PAHFIT (jSmith et al. 
can estimate a total PAH flux (see Table [3]) and compare this to the far-infrared flux 
which we estimate originates from source C+D. When this is done, we derive a ratio of 
FpAHtotl 'Fir ~ 0.02. Since the continuum emission at the overlap wavelengths of SL and 
LL have a difference less than 10%, dust emission in our low-res spectrum is dominated by 
source C+D. 



3.5. AKARI Near-Infrared Spectroscopy 

The spectrum of sources C and D taken with AKARI is well separated from that of the 
other sources (A and B) in II Zw 096, although they are all in the Np aperture, because the 
direction of spectral dispersion is nearly perpendicular to the vector separating C, D, and 
A. The A KARI near- infrared spectrum of source C+D is shown in Figure [3 Strong 3.3 /im 
PAH emission, 3.4 /im emission from aliphatic hydrocarbon grains, and Bra emission lines 
are detected. The line fluxes are (2.37 + 0.14) x 10" 16 Wm~ 2 and (0.50 + 0.05) x lO^Wm" 2 
for the PAH and the Bra, respectively, after fitting both with a Gaussian function (see 
Tables [2] and [3}. The EQW of the 3.3 /im PAH is 0.14 /mi. Because the AKARI spectral 
resolution between 2.5 — 5 /xm is only R ~ 120, the Bra recombination line is not resolved. 

The spectrum of source A (not shown) also exhibits 3.3 /tm PAH and Bra emission. 
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The 3.3 /im PAH has a flux of (1.77 ± 0.12) x lO^Wm- 2 and an EQW of 0.08 /im. The 
Bra line flux is (0.33 ± 0.06) x 10~ 16 Wm~ 2 . 



3.6. Bright Star Clusters In the Remnant 

The high resolution HST imaging data show that a large number of bright star clusters 
surround the II Zw 096 nuclei and are present throughout the merging disks (see Figure |2]). 
From the B- and /-band images, we identify 128 clusters with apparent magnitudes of 
18.3 mag < mg < 26.7 mag (Vega) and 17.9 mag < mi < 25.7 mag (Vega) with 80% 
completeness limits of 26.0 mag and 25.5 mag, respectively. Of these, 97 and 88 are 
detected in the FUV and if-band images, respectively. For all clusters, we use a 0.1" radius 
circular aperture for the HST photometric measurements. We estimate that the clusters 
contribute ~ 14% of the total FUV luminosity (k 8 x 10 36 W) in II Zw 096 (uncorrected 
for reddening). 

In Figure El we present an FUV-optical (F140LP-F435W vs. F435W-F814W) 
color-color diagram of the clusters. We additionally show the colors of a number of diffuse 
emission regions throughout the merger system. For the diffuse emission regions, we use 
0.5" — 0.8" radius circular apertures, depending on the location in the FUV image to avoid 
crowded regions. 

Two potential explanations for the observed color distribution are considered. First, 
if all the clusters are coeval, then their different colors result from variable extinction. 
In other words, those clusters or diffuse emission regions which have colors of F435W- 
F814W > 1.0 mag or F140LP-F435W > 0.0 mag are h eavily reddened (Ay > 3 mag). 



When we de- redden an instantaneous starburst model of iBruzual &: Chariot 



fl2003f ). all of 



the clusters and the diffuse regions are younger than ~ 5 Myr. On the other hand, when 
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we pick an exponentially declining (SF R=lM^r e 



or a continuous star formation model (IBruzual fc Chariot 



' or r = 1 Gyr, < t < 20 Gyr) 



20031 ). they have ages between 



approximately 1 — 100 Myr. However, we note that the points of cluster colors are not 
aligned well with the extinction vector. In particular, there is a significant spread in B — I 
(F435W-F814W) colors, which is not consistent with a coeval population with similar 
intrinsic colors and variable reddening. 

A more likely explanation is that the colors of the clusters imply at least two populations: 
one with F435W-F814W < 0.5 mag and a redder group with F435W-F814W > 1.0 mag. 



The extinction-corrected ages of th e two systems are 



the instantaneous starburst model (IBruzual fc Chariot 



- 5 Myr and 20 - 200 Myr for 



20031 1 . respectively. We assume 



that both have colors consistent with reddening by 1 — 2 magnitudes of visual extinction 
(Ay = 1 — 2 mag). However, many clusters have 0.5 mag < F435W-F814W < 1.0 mag, 
implying ages of 200 — 500 Myr old, if they are un-reddened. If these clusters sit behind 
1 — 2 mag of visual extinction, then they could be as young as 5 Myr old. If we take into 
account an exponentially declining (SF R = lM Pl T~ 1 e^~ t ^ T ' > for r = 1 Gyr, < t < 20 Gyr), 



or a continuous star formation history ( IBruzual fc Chariot 



20031 ) . then the ages of the 
clusters will be between approximately 1 Myr and 1 Gyr with Ay = 1 — 3 mag. With these 
models, there are a larger number of older clusters than in the instantaneous starburst 
model. In particular the clusters that have colors of both of F435W-F814W > 1.0 mag and 
F140LP-F814W ~ -1.0 mag, are about 10 times older. 

In the A ntennae Galaxies , the merger with the most extensively studied system of 
star clusters 



Whitmore et al. 



(119991 ) find three populations of star clusters with ages of 
~10 Myr, ~100 Myr, and ~500 Myr. In the overlap region, the clusters appear to be even 
younger, possibly less than 5 Myr old based on Ha imaging. The youngest clusters in 



II Zw 096 therefore, appear to have similar ages as the youngest clusters in the Antennae. 



- 19 - 



In mo re luminous galaxies (LIRGs and ULIRGs), ISurace et al.l (119981 ) and IScoville et al. 



( 120001 ) find cluster ages between 5 — 300 My r, consistent with the spread in ages we find for 



II Zw 096. In Arp 220, 



Wilson et al. 



(120061 ) similarly find two populations of clusters with 



ages of < 10 Myr and 70 - 500 Myr. 

The regions of diffuse emission in II Zw 096 are typically 0.5 — 1 mag redder in 
F435W-F814W than the star clusters. If we assume that the emission of the diffuse regions 
comes from a mixture of young stars (~ 1 Myr) and old stars (~ 1 Gyr), then the average 
F140LP-F435W color of the diffuse regions can be explained by saying that they contain a 
small fraction (~ 5 — 10%) of young stars. On the other hand, the F435W-F814W color 
of the diffuse emissions suggests that nearly 17% of the emission comes from young stars, 
although the diffuse emissions have a large spread in this color. 

We show the location of the clusters and the diffuse emission on the ACS F140LP 
image in Figure [9j There is no apparent correlation of the cluster colors (ages) with 
position, except for the bluest clusters which are prominent along the spiral arm to the 
northwest. On the other hand, the diffuse emission regions with F435W-F814W < 1.4 mag 
are concentrated in the east while those with F435W-F814W > 1.4 mag are in the west of 
II Zw 096. 

The F435W-F160W or F814W-F160W colors of the clusters do not show an obvious 
gradient or pattern with position. However, the clusters to the east side have less blue 
F814W-F160W colors on average than the clusters to the west. 

The NICMOS if -band data can be used to compute the masses of the clusters. In 
Figure [101 we show an F435W-F814W vs. absolute if -band color-magnitude diagram for 
the clusters. Assuming the clusters have — 2 mag of visual extinction, as in the color-color 
diagram (see Figure |8]), the if-band data suggest that most of the clusters have masses of 
10 6 - 1O 8 M . 



-20 - 



If we compare the masses of the star clusters in II Zw 096 with those in the Antennae 



19991 ). However, we 



Galaxies, the former appear to be ~ 100 times larger (jZhang Sz Fall! 
note that the resolution of the HST data is about 30 pc for II Zw 096 as opposed to only 
5 pc at the distance of the Antennae, and this might account for the apparent differences 
in the cluster mas ses between the two gala xies. For comparison , in (U)LIRGs, the cluster 



masses derived by ISurace et al. 



(|l998h and 



Scoville et al. 



(120001 ) lie between 10 5 - 10 7 M G 



0> 



with the youngest clusters dominating the high mass end. This is also true in II Zw 096, 
where the youngest clusters (those with F435W-F814W < 0.5 mag and Ay = 1 — 2 mag) 
appear to have masses above 1O 7 M . 

In order to compare the distribution of the luminosity of the clusters, we derive 
the .B-band luminosity function (LF) and fit this with a power-law, dcp/dLs oc L B . In 
the range of — 14.7mag < Mb < — 10.7mag with the detection limit of —9.9 mag, we 
obtain a = —1.1 ± 0.3. To avoid confusion due to limited spatial resolution, we chose for 
comparison merging infrared luminous galaxies at similar distances to II Zw 096. The LF 
of the clusters in II Zw 096 is found to have a flatter slope (a = — 1.1 ± 0.3) than that 
of Mrk 231 (a = -1.80 ± 0.2) at z = 0.0 42 or Mrk 463 (a = -1.63 ± 0.11) at z = 0.050, 



which are two nearby ULIRGs studied by 



Surace et al. 



(119981 ) . This suggests that there are 



more luminous clusters in II Zw 096 than Mrk 231 and Mrk 463. In the Antennae Galaxies 
at z — 0.006, the \^-band LF of the clusters have much steeper slopes, a = —2.6 ± 0.2 
and a = —1.7 ± 0.2 in the ranges o f My < —10.4 mag and — 10.4mag < My < —8.0 mag, 



respectively ( iWhitmore et al. 



19991 ). 



4. Discussion 



We have presented the first mid- and far-infrared images and spectra of II Zw 096, whic h 
show that the extremely red source first identified in the near-infrared by iGoldader et al. 
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( 119971 ). is a heavily obscured starburst that dominates the energetics of the entire merging 
system. 



4.1. The Nature of Source C+D 



The IRS spectra of source C+D show the presence of PAH emission, and no evidence 
of [Ne V] or [O IV]. The upper limits on the [Ne V]/[Ne II] and [O IV] /[Ne II] emission line 
flux ratios are 0.04 and 0.11, respectively. These line ratios imply a starburst dominated 
system. The EQW of the 6. 2 jim. PAH feature i s 0.26 /im, which is only about 1/2 that 



seen in pure starburst nuclei (IBrandl et al. 



20061 ) . This might indicate an excess of hot dust 



in source C+D compared to the nuclei of most starburst galaxies, or a deficiency of PAH 
emitters. 

If we adopt a total far-infrared flux Fir = 8.48 x 10~ 13 W m~ 2 for source D, and a 
total PAH flux F PAHtot = 1.44 x 10~ 14 W m" 2 including all the PAH features (Table ED, we 
derive a ratio of F P Am„ t / F rR ~ 0.02. This is at t he low end for pure starburst galaxies, or 
about 1/2 the median value ([Marshall et al.l 120071 ). 



The high [Ne 



starburst galaxies (IThornley et al 



II] /[Ne II] line flux ratio of ~ 1.0 is at t 



2000 



Madden et ai- 



r e high-end o 



2006 



Beirao et al. 



that seen in 



20081 ). There 



are two possibilities that can explain this high ratio - either a low metallicity, or an 
unusually ha rd radiation field, p ossibly produced by a large number of high mass stars. 



According to 



Rupke et al. 



(120081 ). the SE source (source A in this article) has a metallicity 
of 0AZ Q < Z < 1.8Z Q , or a mean oxygen abundance (12 + /og(0/H)) = 8.6. If we assume 
solar metallicity for the gas, the implied age is about 3 Myr, and an upper mass cut off o f 



50 — 100 M Q is required to match the observed [Ne III]/[Ne II] ratio (IThornley et al 



200(1 ) . 



The strong 3.3 /im PAH emission feature seen in the AKARI spectrum also supports 
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a powerful starburst in source D ( IMoorwoodl Il986l ; llmanishi et al.l 120081 ) . In addition , the 



typical indicators o 



Veilleux et al. 



a buried AGN, e.g. a rising (power law) continuum (lArmus et al 



20091 ) or 3.4 fim absorption by aliphatic hydrocarbon grains (llmanishi et al 



2007; 



20061 ) are not detected. 



The Chandra data for source C+D is also consistent with a buried starburst, although 
it is slightly weak in X-rays for its far-infrared luminosit y. Starburst galaxi es show a 



correlation between their infrared and X-ray luminosities ( jRanalh' et al 



20031 ). To convert 



from 2 — 7 to 2 — 10 keV, we ha v e extr apolated the spectral model fit to the data from 



7 to 10 keV as in 



Iwasawa et al. 



( 120091 ). Together, source C+D lie below the starburst 



correlation, but within the 2-cr envelope of the distribution. In some (U)LIRG nuclei of 



the GOALS sample, the low ratio of the hard X-ray e mission to the 



ar-infrared emission 



20091 ). although this 



has been taken as evidence for highly obscured AGN (llwasawa et al. 
argument hinges on the detection of an Fe K emission line. Given the low number of counts 
in the spectrum of source C+D, an accurate description of the X-Ray properties will require 
significantly deeper integrations. The othe r X-ray lumin o us so urce, source A, lies on the 



Ranalli et al. 



(120031 ) . suggesting this source is 



X-ray selected starburst correlation line of 
powered by starburst. 

From the infrared luminosity of Ltr = 6.87 x 1 11 L , the star formation rate of source 



D is estimated to be 120 yr 1 ( 



above that estimated by 



Kennicutt 



Goldader et al 



19981 ). This is approximately a factor of five 



(119971 ) based on the near-infrared data alone. 



Assuming this far-infrared emission originates from a region having a size comparable to that 
of the projected FWHM of the MIPS 24 /im beam (6"), we estimate a starburst luminosity 
density of ~ 4.0 x 10 10 L Q kpc -2 , or a star formation rate density of ~ 6.9 M Q yr -1 kpc~ 2 . If 
instead, we take the size of source D estimated from the NICMOS image (~ 220 pc radius), 
we obtain values which are about two orders of magnitude larger, namely a luminosity 
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density of ~ 4.5 x 10 12 L©kpc~ 2 and a SFR density of ~ 780 M yr _1 kpc~ 2 . Although 
large, this va l ue is s till below the peak values estimated for clusters in starburst galaxies 



Meurer et al. 



(119971 ) which are an order of magnitude higher. Similarly, our derived valu e 



for so urce D is well below most of the ULIRG starburst nuclei observed by 



Soifer et al. 



(120001), except for IRAS 17208-0014. 



The Bra line detected with AKARI provides a direct measure of the star formation 
rate. We assume Case B for hydrogen recombination, effective temperature of ~ 10 4 K, 
density of ~ 100 cm -3 . We then estimate a star- formation rate of approximately 45 M yr _1 



uncorrected for extinction (IKennicutt 



19981 ). Although we do not have an estimate of the 



extinction from the AKARI spectra, if we apply the line of sight value e s timat ed from the 



Smith et al. 



(2007^ eq.(4)). the 



depth of the silicate feature in the IRS data (Ay = 19 mag; 
implied SFR would be about 7OM yr~ 1 . 

Existing radio observation s of II Zw 096 also su ggest that source D is a site of intense 
star formation. At 1.425 GHz, 



Condon et al. 



(119961 ) measure an integrated flux density of 



26.6 mJy using the Very Large Array (beamsize of 6"). The J2000 position of the radio 
peak is R.A.=20 h 57 m 24 s .3, Dec=+17°07'38 8", well within the positional uncertainties of 
the position of source D in the MIPS data. With the n ew Spitzer imagin g data, we can, 



for the first time, estimate the infrared to radio ratio q ( IHelou et al 



19851 ) for source D as 



q = 2.93. This value is withi n the ~ 1-a (0.26 dex) scatter of the average q ((q) = 2.64) 



among star forming galaxies (jBell 



20031 ). and is thus consistent with our mid-infrared 



spectra and X-ray data implying starburst activity in source D. 

The if -band magnitude (—19.5 mag) and assumed ag e, extinction, and metallic ity 



are used to der i ve th e mass of source D by the model of 



Bruzual fc Chariot 



Goldader et al. 



J2003|). 



( 119971 ) estimate the age of the II Zw 096 starburst to be roughly 5 — 7 My r 



from the CO index and Br7 equivalent width, using the models of 



Leitherer fc Heckman 
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( 119951 ). Assuming an extinction of 3.3 mag in if-band (Ay > 19 mag - see section |3~3|) and 



solar metallicity, we estimate the mass of source D to be approximately 1 — 4 x 10 9 M< 



4.2. The Comparison with Other Merger-Induced Extranuclear Starbursts 



Source D is a powerful starburst not associated with the primary nuclei of the II Zw 096 
system (sources A and B). It could be a starburst in the disturbed disk of source A, or 
even the nucleus of a third galaxy - one which is nearly completely hidden behind A and B. 
While large scale mapping and modeling of the gas and stellar dynamics of II Zw 096 will be 



needed to fully understand the true nature of the sys tem, we 



avor the former explanation 



Barnes! ( 12004] )) that shock-induced 



for two reasons. First, simulations have shown (e.g. 
starbursts can occur both in and outside the nuclei of merging galaxies. Second, we have 
a number of examples of powerful, extranuclear starbursts in the local Universe. Below 
we compare the properties of source D to the extra-nuclear starbursts in two well-studied 
merging galaxies NGC 4038/9 (the Antennae Galaxies) and Arp 299. 

In the Antennae Galaxies, nearly 1/2 of the total 15 /im emission (about 5 x 10 10 L Q ) 
comes from the overlap region (5x3 kpc ) between the merg ing spiral galaxies, with about 



15% coming from a single starburst knot (IMirabel et al. 



19981 ). Recent mid- infrared spectral 



mapping observations with the IRS on Spitzer have resolved this off-n uclear starburst reg ion 



into a number of clumps in the mid- infrared (peaks 1, 2, 3, and 5 of 



Brandl et al 



(|2009j)). 



These clumps (or clusters) have infrared luminosities of 3.86 x 10 L@ < Lm < 1.14 x 10 L , 



corres ponding to ~ 5 — 16 % of infrared radiation coming from the system. 



Brandl et al. 



( 120091 ) estimate that the SFRs in these clusters range between 0.66 — 1.97 M yr 1 . 



To make a direct com parison with 



in the Antennae (peak 1 in 



Brandl et al 



I Zw 096, we pick the most luminous infrared peak 
( 120091 )) which has L IR — 1.14 x 1O 1O L and assume 
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the ra dius of this cluster to be R ~ 48 pc (the size estimated by IWhitmore fe Schweizer 
(I1995I )). We then estimate a luminosity density of 1.6 x 1O 12 L kpc -2 . This value is about 
three times smaller than the upper limit of the luminosity density of source D (assuming a 
size of 220 pc radius - the size of D in the NICMOS image). While the luminosity density 
of peak 1 is extremely large, it only contributes ~ 15% of the total infrared emission of 
the entire system, in contrast to source D in II Zw 096 which dominates the total infrared 
output of the system. 

In the A ntennae Galaxies, t he silicate optical depths at 9.7/im of the star clusters 



measured by 



Brandl et al. 



(120091 ) are r 9 . 7 = 0.13 - 1.03 (A v = 2.4 - 19 mag). Although 



one of the clusters in the Antennae, peak 5, has an extinction comparable to that measured 
in source D in II Zw 096 (r 9 7 ~ 1.03), its Ljr is the smallest in all of the clusters in the 
overlap region of the Antennae. This is unlike source D which has both high extinction and 
high Lir. Peak 1 itself, has a r 9 7 = 0.19, much lower than our estimated value for source 



D, even though the estim ated age of pea 



£ 1 (2 — 4 Myr) is slightly less than that derived 



for source D (5 — 7 Myr - iGoldader et al. 



(jl997h ). 



Arp 299 (IC 694 + NGC 3690) has substantial extranuclear infrared emission like 
the Antennae but an infrared luminosity of ~ 6 x 10 11 L Q , making this system a LIRG. 
In particular, regions C and C' in the overlap region betw een IC 694 and NGC 3690 



have all the signatures of a merger-induced starb u rst (e .g. 



Charmandaris et al 



(12002 



luminosity (4.4 x 10 10 L o ; 



Alonso-Herrero et a 



Charmandaris et al. 



Alonso-Herrero et al. 



(120001 ): 



( 20091 )). Region C+C have an infrared 



(120021 ) ). approximately 6.4% of the L IR of 



source D in II Zw 096. The fractio n of the infrared emission in regions C and C to that in 



the Arp 299 system is only ~ 10% (j Charmandaris et al. 



20021 ) , much less than the fractional 



infrared luminosity in source D (80%). The infrared luminosity densities of regions C and 
C are estimated to be ~ 1.0 x 10 11 L Q kpc -2 and ~ 7.5 x 10 10 L kpc -2 , respectively, if 



-26 - 



we assume that their sizes are R ~ 300 pc and R ~ 2 50 pc, corresponding to the size of 



Paa emission estimated by 



Alonso-Herrero et al. 



(120091 ) . and that the infrared luminosity 



of regions C and C have a ratio of 2:1. These numbers are a factor of 2 — 3 larger than 
what we estimate for source D when we assume that the size of the source D is equivalent 
to the beamsize of the MIPS 24 /im. However, if the size of source D is comparable to the 
NICMOS resolution, then the luminosity density could be significantly larger than that 
seen in regions C and C in Arp 299. 



Assuming a metallicity of 2Z Q , 
regions C and C and the models of 



Alonso-Herrero et a 



Snijders et al 



(120091 ) use the IRS spectra of 



( 120071 ) to estimate ages of 4 — 7 Myr. 



If we apply the same model and metallicity to source D in II Zw 096, the measured 
[Ne IIl]/[Ne II] = 1.0 and [S IV] /[S III] = 0.33 line ratios suggest that source D also has an 
age of 4 — 7 Myr. 



5. Summary and Conclusions 

The first analysis of the interacting galaxy II Zw 96 utilizing mid- and far-infrared 
images and mid-infrared spectroscopy leads to the following conclusions: 

1. The Spitzer imaging reveals that approximately 80% of the far- infrared emission of 
the II Zw 096 comes from a compact, off-nuclear starburst with R.A.=20 h 57 m 24 s .34, 
Dec=+17°07'39 1" (J2000). The estimated 8 - 1000 /xm luminosity of this source 

is L IR = 6.87 x 10 11 L Q . The implied star formation rate of this object is about 
^OMoyr- 1 . 

2. HST NICMOS observations show that the off-nuclear starburst is composed of two 
prominent, red knots, with a number of smaller peaks spread over approximately 10 
square arcseconds (~ 6.6 kpc 2 ). Most of this emission is not seen in the optical imaging 
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data, either from the ground or with the HST ACS. The knot to the southwest of 
this region, source D, is consistent with the peak in the Spitzer MIPS 24 and 70 (ttrn 
emission which dominates the far-infrared luminosity of the entire merging galaxy. 
Assuming a standard mass-to-light ratio for galaxies, the H-band luminosity of source 
D corresponds to a mass of 1 — 4 x 10 9 M . 

3. The Chandra X-ray imaging shows that the hardest source in the system is source 
C+D, and it lies slightly below the correlation of X-ray and infrared luminosities seen 
for starburst galaxies, although within the 2-a scatter of the relation. 

4. The Spitzer IRS and AKARI spectra suggest a starburst-dominated system with 
strong 3.3/im and 6.2/im PAH emission, and no evidence of emission from [Ne V] 
14.3//m or [O IV] 25.9/zm in the IRS data. While the equivalent width of the 3.3/zm 
PAH feature as measured in the AKARI spectrum is typical for starburst galaxies, 
the equivalent width of the 6.2/zm PAH feature as measured in the IRS spectrum is 
about 1/2 the value found in pure starbursts. This may reflect an excess of hot dust 
surrounding source C+D, which is easier to detect in the narrow IRS slit. This is 
consistent with the [Ne III]/[Ne II] line flux ratio of ~ 1 indicating a relatively hard 
radiation field. The 9.7/zm silicate optical depth, measured from the IRS spectra 
suggests an A v > 19 mag toward source C+D, indicating a highly buried source. 

5. The ACS B and /-band imaging reveals a large number of (super) star clusters in 
the II Zw 096 system, consisting of at least two populations: one 1 — 5 Myr and one 
20 — 500 Myr, which may have formed at different times during the merger. We find 
no clear association of cluster age with position in the merger. Most of the cluster 
masses are in the range of 10 6 — 10 8 M . 

Spitzer imaging and spectroscopy of II Zw 096 has revealed a powerful, young 
extranuclear starburst. This starburst is reminiscent of those seen in NGC 4038/9 (the 
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Antennae Galaxies) and Arp 299, but it is more luminous (more than an order of magnitude 
more luminous than the Antennae starburst), and it is responsible for nearly all the infrared 
luminosity in II Zw 096 (compared to only 10 — 15% in the Antennae or Arp 299). Source 
D in II Zw 096 is one of the most extreme buried extra-nuclear starbursts yet discovered in 
the local Universe. 
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Fig. 1. — The Spitzer IRAC and MIPS images of II Zw 096. From upper left to lower 
right, these are (a) 3.6 /mi, (b) 4.5 /im, (c) 5.8 /mi, (d) 8 /mi, (e) 24 /mi and (f) 70 /im. In all 



frames, north is up and east is to the left. Sources A, B, C and D in I 



in the 3.6 /im image, using the same nomenclature as 



Goldader et al. 



Zw 096 are labeled 



(119971 ). The origin 

of both axes in each panel corresponds to the position of source D at R.A.=20 h 57 m 24 s .34, 
Dec=+17°07'39 1" (J2000). The contour levels of the 3.6 /mi, 4.5/mi, 5.8 /an, 8.0 /mi, 24 /mi 
and 70 /im images start at 6.0, 9.0, 24, 55, 180, and 100 MJy sr _1 , respectively, and increment 
with a step factor of two. The four IRS slits are overlayed on the IRAC 8 /im image in (d). 
The blue solid, red dashed, orange dotted and green dashed dotted boxes indicate Short-Low, 
Short-High, Long-High and Long-Low, respectively (See the electronic edition of the Journal 
for a color version of this figure). Spectra extracted from the Short-Low and Short-High slits 
are from the peak emission only (source C+D), not from the entire slit, as they are for the 
Long-Low and Long-High slits - see text for details. 
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Fig. 2.— The HST/ ACS and SBC F140LP, F435W, F814W, and NICMOS F160W grey- 
scale images of II Zw 096 with Spitzer MIPS 24 micron contours overlaid. The source 
of the prodigious far-infrared emission is mostly invisible in the optical and UV. The la- 
bels AR.A. = and ADec = indicate the coordinate of source D (R.A.=20 h 57 m 24 s .34, 
Dec=+17°07 / 39 1"). 
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Fig. 3. — Close-up of the NIC MOS image, identifying the red sources designated as "C" and 
"D" by lGoldader et al.l ( 119971 ) . The new HST data, reveal not only a large number of young, 
blue clusters, but also the true, complex nature of the buried starburst. 
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Fig. 4. — False-color RGB image of II Zw 096 made by combining the HST ACS SBC 
F140LP (blue), ACS F435W (green) and NICMOS F160W (red) images. 
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Fig. 5.— The Chandra (a) full band (0.5 - 7 keV), (b) soft X-ray (0.5 - 2 keV), and (c) 
hard X-ray (2 — 7 keV) images. The contour levels of all of the images start at 0.3 with 
a step of 0.2. The labels AR.A. = and ADec = indicate the coordinate of source D 
(R.A.=20 h 57 m 24 s .34, Dec=+17°07'39 1"). 
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Fig. 6. — Spitzer IRS spectra of source D (see Fig [3]). The full low- resolution spectrum is 
at the top left, while a close-up of the Short-Low spectrum is shown at the bottom left. 
The Short-High and Long-High spectrum are shown at the top right and the bottom right, 
respectively. Some residual order tilting is evident in the two reddest orders of Long-High. 
All spectra are shown in the rest frame. 
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Fig. 7. — AKARI near-infrared spectrum of sources C and D. The wavelength scale is given 
in the rest-frame. 
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Fig. 8. — Optical-UV color-color diagram of the young star clusters (black diamonds) and 
diffuse emission regions (red triangles) in II Zw 096. Only clusters whose uncertainties are 
F140LP-F435W < 0.2 mag and F435W-F814W < 0.1 mag are plotted. The maximum uncer- 
tainty is shown at the lower right as an error bar. The solid purple, the orange dotted and the 
blue dashed lines are the models that correspond to instantaneous starburst, exponentially 
declining (SFR=lM Q r~ 1 e^*/' r ^ for r = 1 Gyr, < t < 20 Gyr) and con tinuous star forma- 
tion age tracks, respectively, with Salpeter IMF and Padova 1994 track (jBruzual &: Chariot 
* 0031). The ages are labeled in years. The arrow indicates the Ay = 1 mag reddening vector 
( ICalzetti et al.lll994j ). See Figure [9] to identify the positions of the clusters and the diffuse 



emission regions. 
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Fig. 9.— The HST ACS SBC F140LP (far-UV) image with the positions of the clusters 
and diffuse emissions (see section |3H]) . Blue crosses are the clusters whose F435W-F814W < 
0.5 mag implying the ages are 1 — 5 Myr, red triangles are those which have F435W-F814W > 
1.0 mag implying the ages are 20 — 200 Myr, and green diamonds have colors between 
0.5 mag < F435W-F814W < 1.0 mag. The purple and orange circles indicate the diffuse 
emission and correspond to F435W-F814W < 1.4 mag and F435W-F814W > 1.4 mag, 
respectively. The size of the circles corresponds to the aperture size of the photometry 
measurements. Also see Figure |8j 
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Fig. 10.— Optical-Near-infrared color-magnitude diagram of the young star clusters in 
II Zw 096. Only the clusters which have uncertainties of F160W < 0.4 mag and F435W- 
F814W < 0.1 mag are plotted in the figure. The error bar at the lower right indicates the 
maximum uncertainty. The dotted, dashed and solid lines are the instantaneous star burst 
aging tracks for clusters with 10 8 , 10 7 and 10 6 M Q , respectively. The arrow is the Ay = 1 mag 
reddening vector. 
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Table 1: Spitzer, Hubble, Chandra, and AKARI observations of II Zw 096 



osei vatoi y 


\ Y~i nfv>n TV"! ATI T" 


r lrcei oi mo Quit; 


iircegiaTjion time [secj 


LiDseivation oate 






3.6 /im 


150 






IRAC 


4.5 /im 
5.8 /im 
8.0 /im 


150 
lou 

D 


2004 Oct 29 






24 /im 


48.2 




Spitzer 


1\/TTPQ 


70 /im 
160 /im 


Q7 7 
O ( . ( 


ZUU4 1NOV ou 








U A 14 






IRS 


SH 
LH 


A x 14 a 
6 x 30 a 
4 x 60 a 


2006 Nov 12 




ACS/SBC 


F140LP 


2528 


2008 May 1 


HST 


ACS/WFC 


F435W 
F814W 


2520 
1440 


2006 Apr 15 




NICM0S/NIC2 


F160W 


2304 


2008 Jun 1 


Chandra 


ACIS 


0.5 - 7 keV 


14.56 x 10 3 


2007 Sep 10 


AKARI 


IRC 


NG with Np 


799 


2009 May 12 & 13 



"number of cycles x integration time in each nod position 
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Table 2: Fine structure and molecular lines 



Line 


wavelength [/an] 


flux [xlO- 16 Wm" 2 ] 


Bra 


4.050 


0.50 ±0.05 


H 2 S(3) 


9.660 


0.45 ± 0.05 


[S IV] 


10.511 


0.72 ±0.04 


H 2 S(2) 


12.279 


0.29 ±0.03 


[Ne II] 


12.814 


2.74 ±0.47 


[Ne V] 


14.322 


< 0.10 a 


[Ne III] 


15.555 


2.79 ±0.03 


H 2 5(1) 


17.035 


0.73 ±0.05 


[S III] 


18.713 


2.16 ±0.22 (SH) 


[S III] 


18.713 


3.27 ±0.04 (LH) 


[0 IV] 


25.890 


< 0.29 a 


[Fe II] 


25.988 


0.56 ±0.22 


[S III] 


33.481 


4.77 ±0.07 


[Si II] 


34.815 


4.23 ±0.12 



a 3-a upper limit 



Note. - 



The LH spectra cover not only source C+D but also source A (see text for the detail). 
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Table 3: PAH features 



wavelength [/xm] 


flux [xl(r 16 Wm- 2 ] 


EQW H 


3.3 a 


2.37 ±0.14 


0.14 


6.2 b 


5.57 ±0.19 


0.26 


11.2 b 


3.24 ±0.32 


0.14 


6.2 


14.11 ±0.35 


1.30 


7.7 complex 


83.57 ± 1.14 


7.61 


8.3 


4.55 ± 0.44 


0.35 


8.6 


12.86 ±0.37 


0.99 


11.3 complex 


4.98 ±0.15 


0.16 


12.6 


3.34 ±0.27 


0.05 


13.6 


0.50 ±0.24 


0.01 


14.2 


7.08 ±0.26 


0.09 


17 complex 


13.50 ±0.61 


0.15 



a The AKARI detection. The others are detected by Spitzer. The AKARI spectrum is a mixture of sources C 
and D but the Spitzer spectra represent source D. 

^Measured using a simple spline fit to the continuum, and integrating the feature flux above the continuum. 
Other values are measured using the PAHFIT routine (Smith et al. 2007), which fits Drude profiles to all 
the PAH features, resulting in systematically larger fluxes and EQW 



Table 4. Spitzer, HST, and Chandra Photometries of II Zw 096 



Observatory 


Instrument 


Filter or module 


Source A 
20 h 57 m 24 s .09, +17°07'35 3" 


Source C 
20 h 57 m 24 s .47, +17°07'39 9" 


Source D 
20 h 57 m 24 s .34, +17°07'39 1" 






3.6 fim 


3.29 ± 0.10 (2.9") 


2.26 ±0.08 (1.4") 


2.85 ±0.08 (1.4") 




IRAC 


4.5 fim 
5.8 fim 


2.68 ± 0.08 (2.9") 
8.43 ± 0.18 (2.9") 


3.33 ±0.08 (1.4") 
9.36 ±0.22 (1.4") 


4.67 ±0.13 (1.4") 
12.7 ±0.2 (1.4") 


Spitzer 




8.0 fim 


28.1 ±0.5 (2.9") 


27.2 ±0.4 (1.4") 


41.9 ±0.8 (1.4") 






24 fim 






1.37 ±0.18 a 




MIPS 


70 fim 
160 jim 






10.2 ± 1.2 a 




ACS/SBC 


F140LP 


(6.59 ±0.30) x 10~ 15 (2.9") 






HST 


ACS/WFC 


F435W 


(3.03 ±0.02) x 10~ 15 (2.9") 




(3.43 ±0.34) x 10~ 17 (1.4") 
(4.63 ±0.67) x 10" 18 (0.25") 




F814W 


(1.28 ±0.01) x 10~ 15 (2.9") 


(5.57 ±0.19) x 10~ 17 (1.4") 
(4.64 ±0.39) x 10" 18 (0.25") 


(4.85 ±0.20) x 10~ 17 (1.4") 
(2.85 ±0.08) x 10" 17 (0.25") 




NICMOS/NIC2 


F160W 


5.36 ± 0.05 (2.9") 


0.85 ±0.02 (1.4") 
0.13 ±0.01 (0.25") 


0.74 ±0.02 (1.4") 
0.24 ±0.01 (0.25") 


Chandra 


ACIS 


0.5 - 2 keV 
2-7 keV 


6.6 x 10 40 
1 x 10 41 


1.8 x 10 40 (Source C+D) 
1 X 10 41 (Source C+D) 



Note. — The values in the brackets are the circular aperture sizes in radii. The units are mjy for IRAC and NICMOS/NIC2, Jy for MIPS, 
ergs" 1 cm" 2 A" 1 for ACS/SBC and ACS/WFC, and ergs" 1 for ACIS. 

a Measured by the PSF fitting (see fl3H . 



